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Abstract. We introduce a Monte Carlo scheme for the approximation of the solutions of fully 
non-linear parabolic non-local PDEs. The method is the generalization of the method proposed 
by [Fahim-Touzi-Warin, 2011] for fully non-linear parabolic PDEs. As an independent result, we 
also introduce a Monte Carlo Quadrature method to approximate the integral with respect to Levy 
measure which may appear inside the scheme. We consider the equations whose non-linearity 
is of the Hamilton-Jacobi-Belman type. We avoid the difficulties of infinite Levy measures by 
truncation of the Levy integral by some k > near 0. The first result provides the convergence 
of the scheme for general parabolic non-linearities. The second result provides bounds on the rate 
of convergence for concave non-linearities from above and below. For both results, it is crucial to 
choose k appropriately dependent on h. 

1. Introduction 

The present paper generalizes the probabilistic numerical method in [10] for approximation of 
the solution of fully non-linear parabolic PDEs to non-local PDEs (integro partial differential 
equations). Here by non-local PDE, we mean the integro-partial differential equations which 
sometimes are simply referred to as IPDEs. The method is originated from [8] where a similar 
probabilistic numerical method is suggested based on 2BSDEs, i.e. second order backward stochastic 
differential equations. 

A probabilistic numerical (Monte Carlo) method is introduced in |10j for fully non-linear par- 
abolic local PDEs. A convergence result is established for a class of non-linearities and bounds 
for the rate of convergence for concave and Lipschitz non-linearities for the Monte Carlo method 
is derived. Moreover, it is shown that the error due to estimation of the conditional expectations 
does not change the results if enough number of sample paths are used. Finally, some numerical 
results are provided for mean curvature flow problem and two and five dimensional continuous-time 
portfolio optimization problems. 

Fully non-linear PDEs arise in many problems in applied mathematics and engineering including 
finance. For example the problem of motion by curvature, portfolio optimization under different 
type of constraints, option pricing under illiquidity cost, etc. Non-local fully non-linear PDEs 
arise from stochastic optimization problems for controlled jump-diffusion processes e.g. problem 
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of portfolio optimization in Levy markets. There are only few examples with explicit and quasi- 
explicit solution; for example [3] and [4]. In some applications the dimension of the PDE is so large 
that classical algorithms like finite element and finite difference fail to approximate the solution 
in a reasonable time. The advantage of Monte Carlo methods is that they are less sensitive with 
respect to dimension on comparison with other methods. 

As in [TO] , the main idea is to separate the equation into a purely linear part and a fully non-linear 
part. Then, we use the time discretization of a suitable jump-diffusion process to approximate the 
derivatives and integral term in the non-linear part. The separation into linear and non-linear part 
is arbitrary up to the satisfaction of some assumptions. The assumptions needed for this result are 
degenerate ellipticity condition for the remaining non-linearity and that the diffusion coefficient is 
needed to dominate the partial gradient of the remaining non-linearity with respect to its Hessian 
component. 

The other contribution of this paper is the Monte Carlo method for approximation of the integral 
with respect to Levy measure which appears in the non-local PDEs. The method is referred to 
in this chapter as Monte Carlo Quadrature (MCQ). We treat the jumps as in [7j for finite activity 
jump-diffusion processes. For infinite activity jump-diffusion processes, we truncate the Levy 
measure near zero and then treat them as in the finite measure case. We introduce bounds for the 
truncation error with respect to the derivatives of integrand and truncation level. 

Although MCQ is independent of the numerical scheme, we choose to approximate the Levy 
integral inside the non-linearity by MCQ. In this case, we also need to choose appropriate truncation 
bound with respect to time step which retains the convergence and rate of convergence as in the 
local case in [TP] , 

The idea of the proof is captured from [JJ and [TO] for the convergence result and from [3] and [TO] 
for the rate of convergence. However, in the non-local PDEs, we need to conquer the new difficulties 
due to lack of Lipschitz continuity of non-linearities appearing in many interesting IPDEs e.g. HJB 
equations. More precisely, if Levy measure inside the non-local integral is not finite, then the 
non-local non-linearity is of HJB type will not be Lipschitz. This difficulty makes it impossible to 
use directly the method in [TO]. We show that if the truncation level; k; is properly dependent on 
time step; h; then one can produce the approximate solution which converges to the solution of the 
non-local problem. 

The first result concerns the convergence of the approximate solution obtained from the Monte 
Carlo scheme to the viscosity solution of the final value problem. As mentioned before, infinite 
Levy measure in the non-linearity makes the direct use of the method in [TO] impossible due to 
the lake of Lipschitz continuity of the non-linearity. Moreover, if we truncate the Levy measure, 
the non-linearity is Lipschitz but as truncation level tends to 0, the Lipschitz constant blows up. 
We solve this problem through manipulation of the original final value problem to another one 
whose corresponding scheme is monotone. Turning the manipulation back, we obtain a bounded 
approximate solution. This approximation is near the approximate function created by the proposed 
scheme, if the truncation level depends appropriately on h. 
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The second result provides the rate of convergence in the case of concave non-linearity. The 
proof of the rate of convergence uses the results in [5j and [6] which generalizes the result of [T] to 
non-local case. The method is based on approximation of the solution of the equation with regular 
sub and super-solutions. Plugging the regular sub or super-solution into scheme and then usage of 
consistency provide the upper and lower bounds. Here, we also need to impose the condition that 
truncation level depend appropriately on the time step in order to preserve the rate of convergence 
after truncation. For the rate of convergence, we also need to manipulate the equation to obtain 
strict monotonicity for the scheme which is a crucial requirement in using the method in [SJ. 

Finally, as mentioned in |10] for non-local case, it is worthy of noticing the relation with the 
generalization of [T2] to non-local case introduced in [TT] which provides a deterministic game 
theoretic interpretation for fully non-linear parabolic problems. The game consists of two players. 
At each time step in a predetermined time horizon, one tries to maximize her gain and the other to 
minimize it by imposing a penalty term to her gain. More precisely, she starts in an initial position 
and chooses a vector p, a matrix T, and a function <p. Then, he will plug an arbitrary vector w 
together with p, T and (p in a non-linear penalty term which she should pay to be allowed to change 
her position by taking one step with appropriate length in the direction of vector w. At the final 
stage, she will earn as much as a function of her final position. As time step goes to zero, her value 
function at any time and any position will converge to the solution of a fully nonlinear parabolic 
IPDE whose non-linearity relates to the penalty term. Vector p, a matrix T and a function ip 
represent the first and second derivatives and the solution function, respectively. 

The parer is organized as follows: In Section [21 the problematic features of non-local fully non- 
linear PDE is discussed on a nallve generalization of the Monte Carlo method from local case in 
[TO] to non-local case. In Section [3] the Monte Carlo quadrature (MCQ) is presented as a purely 
Monte Carlo approximation of with the error analysis. Section U] contains the results of convergence 
and asymptotic properties of the scheme. 

Notations For scalars a, b S K, we write a Ab := min{a, b}, aV6 := max{a, b}, a~ := max{— a, 0}, 
and a + := max{a,0}. 

By M(n, c2), we denote the collection of all n x d matrices with real entries. The collection of all 
symmetric matrices of size d is denoted and its subset of nonnegative symmetric matrices is 
denoted by Sj". For a matrix A £ M(n, d), we denote by its transpose. For M(n, d), we 

denote A - B := Tr[A T i?]. In particular, for d = 1, A and B are vectors of W 1 and A ■ B reduces to 
the Euclidean scalar product. 

We denote by Cd, the space of bounded continuous functions from [0, T] to M. d . For a suitably 
smooth function ip on Qt ■= (0, T] x M, d , we define \ip\oo '■= su P(t,x)£Q T l9 9 (^ 2; )l an d Mi := Moo + 

\<p{t,x)-<p(t',x')\ 

sup r- 

QtxQ t (x-x>)+\t-t'\? 
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2. Preliminaries and features for non-local PDEs 

Let fi, a be functions from [0, T) x R d to R d and M(d, d) and n be a function from [0, T] x R d x E d 
to R d and M(d, d) and a = £T T (7. Suppose the following non-local Cauchy problem: 

-C x v{t, x) — F (t, x, v(t, x),Dv(t, x), D 2 v(t, x),v{t, •)) = 0, on [0, T) x R d , (2.1) 

v(T,-)=g, on G R d . (2.2) 

where F : R + x E d x E x R d x S d x C d ->■ E and £ x given by: 



+ y ^¥>(£, ^ + x, z)) - ip(t, x) - l { | 2 |< 1} L>(/?(t, x)r](t, x, z)j du(z). 

C x is the infinitesimal generator of a jump-diffusion, Xt, satisfying SDE: 

dX t = fi(t,X t )dt + a(t,X t )dW t + / T)(t,X t -,z)J(dt,dz)+ / r}{t,X t -,z)J{dt,dz), 

J{\z\>l} J{\z\<l} 

where J and J are respectively a Poisson jump measure and its compensation who associate to 
Levy measure v by: 



v(A) = E 



J([0,l],dz) 



J(dt,dz) = J(dt,dz) — dt x v{dz). 

For more detailes on jump-diffusion processes, see [2] and the references therein or the classic work 
of [13]. 

The classical solution for the problem (|2.ip ~ (|2.2p does not exist in general and therefore we 
appeal to the notion of viscosity solutions for non-local parabolic PDEs. 

Definition 2.1. • The viscosity sub (super) -solution of \2.1\) - l2Jty is a upper semi-continuous 
(lower semi-continuous) function y_(v): [0,T] x R d — > R such that: 

(1) for any (toi^o) ^ [0^) x an d an V smooth function ip with: 

= max(min)-{j; — ip} = (v — (p)(to, xq) 

We have: 

>(<) -£ x <p(to,x )-F(;<p,D<p,rP<p,<p(-)) (to,x ). 

(2) g(.)>v(T,-)(<v(T,-)). 

The function v which is both viscosity sub and super solution, is called viscosity solution of \2. 1}) - 

• We say that \2. 1\) has comparison for bounded functions if for any bounded upper semi-continuous 
viscosity super-solution v and any bounded lower semi-continuous sub-solution v, satisfying 

v{T, ■) > v(T, •), 

we have v > v on [0, T] x R d . 
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2.1. Discretization of the jump— diffusion process. Our purpose is to introduce a Monte Carlo 
method which approximates the solution of problem (I2.ip -( j2"?2"j) . For this purpose, we first need to 
provide a discretization for the process X. 

Suppose that h = ^, U = ih, and n > 0. We define the Euler discretization of jump-diffusion 
process X t with truncated Levy measure by: 

Xj; x ' K = x + fl(t,x)h + a(t,x)W h + / r]{t,x,z)J([0,h},dz), (2.3) 

J{\z\>k} 

*K = -Vl" V; ' - K and = x. (2.4) 

where jl(t, x) = fi(t, x) + Jn^xy ??(£, x, z)v(dz) and we make the choice of n = when v is a finite 
measure. Let iV* and be respectively the Poisson process derived from jump measure J by 
counting all jumps of size greater than k which happen in time interval [0, t] and its compensation, 
i.e. 



N t K = / J([0,t],dz) and N t K = / J([0,T],dz). (2.5) 

J{\z\>k} J{\z\>k} 

One can write the jump part of X^ X ' K as a compound Poisson process (see for example [9]) 

= x + /i K (t,x)/ l + f j(t,x)W h + ^r ? (t,x,Z i ), (2.6) 

i=i 

where fi K (t,x) = n(t,x) — J*{ K< | Z | <:L } ^(^ z)v(dz), ZjS are i.i.d. R* — valued random variables, 
independent of W and iV K , and distributed as l^ z \ >K yj-u(d. 



z . 



2.2. The scheme for non local fully non linear parabolic PDEs. In this section, we intro- 
duce a probabilistic scheme by following directly the same idea as the scheme for the local PDEs. 
Then, we consider some problems which prevents us to utilize the scheme in many interesting ap- 
plications. Therefore, we introduce a modified version of the scheme which works for the class of 
non-linearities of HJB type (Hamilton-Jacobi-Bellman). 

Following the same idea as in |10j . one can obtain the following immature scheme. 

v h (T,.)=g and v h (U, x) = T h [v h ](ti, x), (2.7) 

where for every function ip : R + x R d — > R with exponential growth: 

T h [^](t,x):=EU(t + h,X t f ' x )] +hF h (t,x,V h ^,iP(t + h,-)), (2.8) 



where 



where 



V h 1> := (V^Vl^Vfo) , 
V k h ijj(t,x) := E[4(t + h,X t * K )H£(t,xj\,k = 0,1,2, (2.9) 

"I W h uh _ ( W h W% ~ hl d _ x 



a 



h ' 2 K > h? 
The details of approximation of derivatives with (j4.5|) can be found in Lemma 2.1 in fin 
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We intend to extend the result of [lOj to the non-local case. First observe that there is an obvious 
extension which could be done immediately by adding the following assumptions to Assumption F 
in [Eg, i.e. 

Assumption F (i) The nonlinearity F is Lipschitz- continuous with respect to (x, r,p, 7, tp) uni- 
formly int, and \F(-, •, 0, 0, 0, 0)|oo < 00; 

(ii) F is elliptic and dominated by the diffusion of the linear operator C x , i.e. 

VyF <a on R d xRxl' , x§ d xC d ; (2.10) 

(hi) F p e Image(F 7 ) and \F p T F~ F p \ oQ < +00. 

We remind that the non-local non-linearity F is called elliptic if 

(1) F is non-decreasing on the second derivative component, i.e. 

F{t,x,r,p,i u i{)) < F(t,x,r,p,-f 2 ,ip) for 71 < 72- 

(2) F is non-decreasing on the non-local component, i.e. 

F(t,x,r,p,j,i/ji) < F(t,x,r,p,j, ip 2 ) for Vi < V>2- 
Then we have the following Theorem. 

Theorem 2.1. Let Assumption F in hold true, and \^\\, |cr|i < 00 and a is invertible. Also assume 
that the fully nonlinear PDE l\2.1\) has comparison for bounded functions. Then for every bounded 
Lipschitz function g, there exists a bounded function v so that 

v h — > v locally uniformly. 

In addition, v is the unique bounded viscosity solution of problem \2. l]) - [2~^) . 

The proof is an straight forward implementation of the Subsection 3.2 of |10| . 

One of the major class of fully non-linear PDEs, is the class of HJB equations which come 
from stochastic control problems arising in many application including finance. However, The non- 
linearity of HJB equations do not satisfies Assumption F in general. Even for local PDEs of HJB 
type, Assumption F is not valid, because F is not uniformly Lipschitz with respect to x. In addition, 
when the Levy measure u is an infinite Levy measure, there is no chance for F to be uniformly 
Lipschitz within respect to ip. Therefore, we need to develop another theory for HJB equations. 

The other problem which occurs in many applications is the lack of explicit form for non-linearity 
F. We present the following example in order to mention this problem. 

Example 2.1. Suppose we want to implement the scheme for the fully non-linear equation of the 
form: 



-v t - F(x, Dv(t, x),D 2 v(t, x),v(t, •)) 



= 
= SO), 



A PROBABILISTIC SCHEME FOR FULLY NON-LINEAR NON-LOCAL PARABOLIC PDES 



7 



where 

F(x,p,7,^) := sup { £ e (p,7) + / ip{x + 6z)v{dz) \ (2.11) 



C e (p, 7 ) := 06p+l0 2 a 2 7 (2.12) 
l(x,ip) e := / ^(x + 02!)i/(cfa). (2.13) 

TTiis /uZZy non-linear equation solves the problem of portfolio management for one asset in the 
Black-Scholes model including jumps in asset price. For the sake of simplicity, for the moment we 
forget about infinite activity jumps. Observe that if v = (the asset price do not jump) then F 
becomes of the form: 

F(x,p,j,ip) := sup \obp+-6 2 a 2 ^ 
which could be given in explicit form by: 

(bp) 2 



2aV 



and the scheme could be easily implemented as in |10j as well as more complicated examples. But, 
when v 7^ (jump do exists), the explicit form for F is not known and the supremum should 
be approximated. This problem is in common with other numerical methods for fully non-linear 
PDEs e.g. finite difference. Although his problem is obviously beyond the subject of this paper, we 
addressed it in this paper in order to mention that why we need to approximate the integral inside 
the supremum. 

The other problem, which appears in high dimensions, is the calculation of Levy integral inside 
supremum. Some numerical methods to approximate the supremum based on the calculation of the 
linear operator C e +I 8 inside the supremum for different 6s. Therefore, we proposed a Monte Carlo 
Quadrature method to approximate the integral in a purely probabilistic way. The MCQ could be 
considered independently in other applications. 

From now on, we relax the assumption that v is a finite measure. To be precise, we need to 
suppose that \2.13) is of the form 

T{x,f) := / (il)(x + 6z)-^{x)-l { \ z \< l} 0Dil)(x)- z)v(dz). 
Jr* 

In this case, there are two ways to work with singular Levy measure in numerical experiments; one 
is to truncate Levy measure near zero (as we also did for discretization of X) and the other is to 
approximate infinite small jumps by a Brownian motion. In both cases, the general form for the 
approximate F is 

F K (x,r,p,j,ip) := sup < c K r + 9b K p + \o 2 a 2 -f + [ ip(x + Qz)v{dz) > . 

2 h\A»i\ 
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where 

c K := / v{dz) and b K := b / zv(dz). 

J{\z\>k} ^{1>|z|>k} 

We will introduce the modified scheme (|4.3j) in Section |4] based on the approximation of non- 
linearity F obtained from truncation of infinite Levy measure and MCQ and then provide asymp- 
totic results as in [10] for non-local case. 

The generalization of the result of [10] for non-local PDEs would be easy if the function F K were 
Lipschitz uniform on k. But, for infinite Levy measures, this is never the case. To overcome this 
problem, we will show that k could be chosen dependent on h, so that the corresponding scheme 
satisfies the requirements of [T] for the proof of convergence. 



3. Monte Carlo Quadrature (MCQ) 
In this section, we propose a Monte Carlo method the value of the following Levy generator: 

I[tp](x) := / (ip(x + r](z))-ip(x)-t { \ z \< 1} ri(z)-Dip(x))i;(dz). (3.1) 



The method is pure Monte Carlo method to approximate (J3JJ) and, therefore could be used in 
the approximation of Levy integral inside the scheme (|4,3p . Because, the result of this section is 
independent of the numerical scheme (|4.3p introduced in this paper, we organize this Section so 
that one can read it independently from other Section. 

Through out this Section, we drop the dependency with respect to (t, x) or other variables and for 
the sake of simplicity and just write r)(z). (For example in assumption F in Section U rj a, ^(t, x, z) 
which depends on (t, x, z, a, 0) will be considered as rj(z)). 

Notice that in order for (|3.ip to be well-defined for regular functions, we impose the following 
assumption on r\: 

\n(z)\ 

v 71 < C, for some constant C. (3.2) 
|z|Al ~ { ' 

We present MCQ in three cases with respect to the behavior of Levy measure near zero: 

• finite measure; J/ui<n v{dz) < oo, 

• infinite measure; 

- case I: f{\ z \<xy \v( z )W(dz) < °°> 

- case II: Jn z \<i\ \v( z )\ 2 v(dz) < oo. 

3.1. Finite Levy Measure. When Levy measure is finite, we choose k = 0. In this case, we 
introduce Lemma 13.11 which proposes a way to approximate the Levy integral of general form: 



<p(x + ri(z))C(z)dv(z), (3.3) 
and then we use this Lemma to approximate the Levy infinitesimal generator (|3,ip . 
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Let J be a jump Poisson measure with intensity given by Levy measure v, and {Nt}t>o be the 
Poisson process given by N t = J * f Rd J(ds,dz) whose intensity is A := f Rd v{dz). By (I2.6p . we can 
write X x by 



X? = x 



(3.4) 



i=l 



where Z{S are i.i.d. random variables with law j-u(dz). We also introduce a Levy process Yt by 



(3.5) 



i=i 



Next Lemma shows that (|3.3p could be approximated by a Monte Carlo formula purely free of 
integration. 



Lemma 3.1. Let 



uf(p){x) :=E 



Then, for every bounded function ip : 



Proof. For the sake of simplicity, we just concentrate on the jump part of process X x and without 
loss of generality, we write X? = x + X^;=i r ?(^)- The right hand side can be expressed as: 



7J(Z))C(*)<M*) 



(3.6) 



E 



-Xh 



oo 



Y h \N h = n 



n=0 



(\hy 



nl 



Then by (|3^H33|) . 



E 



<p(Xh)Y h 



-Xh 



n=l 

oo 



(A/i) 



n-l 



/;.! 



-A/i 



-l ™ 



n=l 



99 x 



8=1 



Notice that in the above expression, the summation starts from n = 1 because 1^ = when TV^, = 0. 
Because Z;s are i.i.d. one can conclude that, 



Then, one can write 



i=i 



nE 



i=l 



E 



V x 



+ V (z 1 )+j2v(z i ))azi 



8=2 



^h(Z) + J^C(^)|JVh = "-l 
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where Z is dependent of ZiS but has the same law az ZjS. Therefore, we can conclude that: 

E[<p(T,(Z)+XZ)«Z)\N h = n-i 

n=l 



E 



e~ Xh \h 



But, we know that 

oo 

Y,®[p(v(Z) + X%)aZ)\N h 



(n- 1)!' 



-Ah 



71—1 



71=1 



(Xh) 



n-l 



(n- 1)! 



E 



<p(r,(Z)+X^C(Z) 



Therefore, 



E 



<p(X%)Y h = \hE{v(rj(Z)+X* h )C(Z)}- 



Because the density of Z is v ^ dz ' 



A ' 



E 



WE 



<p( V (z) + X£)((z)dv(z) 



□ 



In the light of Lemma (|3.ip . we propose the following approximation for (|3. 1 [) : 

X h [(p](x) := u^ 1 - ip(x) \ u(dz) - Dif(x) ■ / rj(z)v(dz). 

Next Lemma provide error bound for this approximation. 

Lemma 3.2. For any Lipschitz function ip we have: 

|(Z ft -Z)[¥>]|oo < cVfc|£M«. (3.7) 

Proof. As a direct consequence of Lemma (|3TTjl . z^' 1 = iE[^(X^)iV^. Therefore, one can conclude 
that, 



So, because 



E 



< C [h / |7/(«)|l/(d«)+V^ , 



(3. 



which provides the result. 



□ 



3.2. Infinite Levy Measure. In the case of singular Levy measure, we truncate Levy measure 
near zero and reduce the problem to a finite measure. In other words, for any k > we have the 
truncation approximation of integral operator (|3. 1 1) . 



l K [tp](x) := (ip(x + r}(z))-(p(x)-l { \ z \< 1} r}(z)-Dip(x))v(dz). 

J{\z\>k} 

Then, we use Lemma (|3.ip to present the MCQ approximation for (|3.ip . 



l Kjh [(p](x) := - <p(x) dv{z)- I rj(t, x, z) ■ Dtp{x)du{z), 

J{\z\>k} J{1>\z\>k} 
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where by Lemma (|3.ip 

vf h := / <p(X^ + V (t, x, z))v{dz) = h^E Ux x h ' R m 

J{\z\>k} L 

Following Lemma provides the error of MCQ approximation of ()3. 1 j) in the case of infinite Levy 
measure. 

Lemma 3.3. Let function ip be Lipschitz. 
(1) 7 //{| z |<i} \zW(dz) <cc, then 



(3.9) 



\(l Kjh -l)[v]\ 00 <C\Dip\ 00 [Vh+ [ \z\u(dz)). 

\ ^{0<|2|<k} J 

(2) // /rui^i \z\ 2 v(dz) < oo, then 

\(T R>h -T)[y]\^<c(\Dip\JVh + hf \z\u(dz)) + \D?<p\ O0 [ \z\ 2 u(dz)). (3.10) 

V V J{\z\>k) ' J{0<\z\<k} J 

Proof. 

(1) Notice that, 

By (|3.2p . the truncation error is given by: 

|(X -X«)M |oo < 2|Z>v»|oo / M*)K<fc)- (3-11) 

J{0<|z|<k} 

On the other hand, by (|3.8f) and (|3.2I) . we observe that 

KXk-X^JMIoo < ClDipl^lhf \ V (z)\v(dz) + Vh) 

< ClDipl^lh I \z\v{dz) + \[h ) 

which together with (|3.1ip provides the result. 

(2) By (|3.2p . the truncation error is given by: 

|(Z-Z«Moo<C|DV|oo /" |z|M^), (3.12) 

J{0<|z|<k} 



for any function with bounded derivatives up to second order. On the other hand, ([3? 
allows us to calculate the Monte Carlo error by: 

1(2* -:£«,*)[¥>] loo < C\Dip\Ah\ \z\u(dz) + Vh) 

V •ni*i>«} / 

which completes the proof. □ 
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4. Asymptotic results 

This section is devoted to the convergence result for the scheme (|4.3h . We first remind the notion 
of viscosity solution and provide the assumptions required for the main results together with the 
statement of main results. Then we provide the proof of the results in two following subsection. 

we need to impose the following assumption on the non-linearity F to obtain the convergence 
Theorem. 



Assumption IHJB1: Function F satisfies: 

-a(t,x) -7 + /x(t,x) ■ p + F(t,x,r,p,j,^):= inf ^sup {c a > p {t, x, r,p, 7) + l a '^(t, x, r,p, 7, 

for given sets A and B where 

£ a >P{t, x, r,p, 7) := ^a a >P(t, x) ■ 7 + b a ^(t, x) ■ p + c a ^(t, x)r + k a ^{t, x), 

and 

l a ^(t, x, r,p, := (x + rf'^t, x, z)) - r - l {lz ^ 1}V a ^ (t, x, z) ■ p) u(dz) 
where for any (a,/3) G A x B, a a,/3 , b a ^ , c a, @ , k a ^ and rj a ^ satisfy 

sup {\a a % + \b a >% + \c a '% + \k a >% + ^fi-i^h 1 < QO _ 
a&A,0£B I \z\ A 1 J 

The non-linearity is dominated by the diffusion of the linear operator C x , i.e. for any t, x, z, a 
and j3 

\a~ ■ a a *^* |i < 00 and < a a 'P < a, (4.1) 
jffi^ \f*£ G Imaged) and sup |(6 Q ' /3 ) T (a a ' /3 )-6 a ' /3 | 0O < 00, (4.2) 

sup j— ; < OO 

sup —75 < OO. 

aeA,peB l/\\z\ z 

Remark 4.1. A function F which satisfies Assumption IHJB1 is not well-defined for arbitrary 
(t, x, r, p, 7, ip) G R_|_ xl^xlxR^xS^x But, for any second order differentiable function, ip, with 
bounded derivatives with respect to x, F(t, x, ip(t, x), Dtp(t, x), D 2 ip(t, x), tp(t, •)) is well-defined. 



Now, we propose a Monte Carlo scheme for (|2.ip - (|2.2p based on the same idea as in [10] . and 
also the approximation of the non-linearity. 

v K > h (T,.)=g and v K >\t h x) = T K:h [v^ h ](t h x), (4.3) 

where for every function ip : R + x M. d — > 1R with exponential growth: 

T^ h ty](t,x):=E\i/,(t + h,Xt*' K )] +hF Kth (t,x,T> h i/j,^(t + h,-)), (4.4) 
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F K>h (t, x, r,p, 7 , V) = inf sup(^(i 5 x) ■ 7 + b a >?(t, x) ■ p + c a ^(t, x)r + k a >?(t, x) 

^ (vt'+m, -)){x) - r - rf^ x, z) ■ p^j v{dz) } , 



/{|*|>K}\ 



+ 

and 

T>fy(t,x) := E^(t + h,X t f ; x > K )H%(t,x)\,k = 0,1,2, (4.5) 

where 

H - 1, H 1 - [a ) — , # 2 - (ct J ^ cr . 

The details of approximation of derivatives with (|4.5h can be found in Lemma 2.1 in [10J. In order 
to have the convergence result, we also need to impose the following assumption over F K> h- 
Assumption Inf— Sup: For any k > 0, t £ [0, T], x and x' G R d and any Lipschitz functions if) 
and ip, there exists a (a*, (3*) £ Ax B such that 

*?>P [if,, v] (t, x, x 1 ) = jf'P* [m, x) - jf'P M (t, x') 

where 

■ a / ty, <p] (t, x) := inf J?f> M(t,x)- sup [ip}(t,x% (4.6) 

a p 



and 

+ 



?'^](t,a;) := ^a a>/3 • D 2 <j)(t,x) + 0^ • D<j)(t, x) + c a '^(t, x) + k a ^{t,x) 

(i>f "'^(t, .))(z) - Ht, x) - n a '?(t, x, z) ■ D<j>(t, x)) v(dz). 



{\z\>k} 



The first result concerns the convergence of the convergence of v K,h for k appropriately chosen 
with respect to h. 

Theorem 4.1 (Convergence). Let 77, fi and a be bounded and Lipschitz continuous on x uniformly 
on t and z, a is invertible and Assumptions IHJB1 and Inf— Sup hold true, and assume that \2. 1\) 
has comparison for bounded functions. Then, if Kh is such that: 

lim Kh = and lim sup 6*j? /i = (4-7) 



where 



with 



sup 1 9 a J I oo, (4i 



I v(dz) + -(b^-f ^{z)v{dz\ 

J{\z\>k} 4 V ^{l>k|>K} 

x (a ' 15 )- (b a 'P- [ r) a 'P(z)u(dz)) , 
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then v Kfl,h converges to some function v locally uniform. In addition, v is the unique viscosity 
solution of fPl)-fP). 

Specially, if Levy measure is finite for the choice of = the assertion of the Theorem hold true. 

Remark 4.2. It is always possible to choose k/j such that (|4,7p is satisfied. To see this, notice that 
9 K in (14. 8h is non-increasing on k 

lim 6 K = +oo and lim sup # K < oo. 

Then, we define Kh := mf{K\8 K < h~z} + h. By the definition of Kh, 6 Kfi < h~* . Because Observe 
that Kh is non-decreasing with respect to h and lim^o = 0. 

If there exists a q such that, q := lim^o K-h > 0, then, for k < q, we would have 6 K = oo which 
obviously contradicts the fact that for k > 0, 8 K < oo. Therefore, satisfies (|4,8p . 

Remark 4.3. The choice of in the above Theorem seems to be crucial for the convergence. 
Otherwise, we only have the following convergence result. 



Proposition 4.1. Under the same assumption as Theorem \4-l\ when Levy measure v is infinite, 
for every Lipschitz bounded function g, we have 

lim lim v K,h = v 

where v is the unique viscosity solution of \2.1\ )- l2Jty assuming that it exists. 

Proof. . Let v K be the solution of the following problem: 

-C x v K (t,x)-F K (t,x,v K (t,x),Dv K (t,x),D 2 v K {t,x),v K {t,-)) = 0, on[0, T)xR d , (4.9) 
v K (T,-)=g(-), on G R d . (4.10) 

where F R : R + x R d x R x R d x S d x C d R is given by: 

F K (t,x,r,p,j,%l>) := inf 'sup \C a 'P{t,x,r,p,'y) + !%>P(t,x,r,p,'y,iJ))} 

where 

X^(t,x,r,p, 7 ,#=/ (V(x + ^(t,x,z)) -r-l {W < 1} ^(*,x,z)-p)i/(^) (4.11) 
J{\z\>k} 

where a a,/3 , b a ^ , c"'' 3 , /c"'' 3 and r/"'' 3 are as in Assumption IHJB1. Let v K,h be the approximate 
solution given by the scheme (|4.3|) . Let k > be fixed. Because the truncated Levy measure is 
finite, by Theorem 14.11 v K,h converges to v K locally uniformly as h — > 0. Let v K be the solution of 
T9|) - (|4.10p . By Theorem 5.1 of [6] and Assumption F, we have: 

k-^U < C7supi(/ \n a 'P(;z)\lv(dz) ] J> (4.12) 

a,/3 \Jo<\z\<K 




< C[ I \z\lv(dz) . (4.13) 

'0<|2|<K / 
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Therefore, one can choose k > so that \v K — v\oo be small enough. Then, when h goes to 0, v K,h 
converges to v K . □ 

The above limit proposes to implement the numerical scheme in two steps: 

• First by choosing k so that v K is near enough to u, we obtain a uniform approximation of v. 

• Second by sending h — > 0, we obtain locally uniform convergence of v K,h to v K . 

Notice that the above convergence is not uniformly on (k, h). However, the convergence in Theorem 
14. 14 is uniform on h when the choice of k is made suitably dependent on h. 



Remark 4.4. By Remark 3.7 in [ID], the boundedness condition on g can be relaxed. 

In order to obtain the rate of convergence result, we impose Assumptions IHJB2 and IHJB2+ 
which restrict us to concave non-linearities. 

Assumption IHJB2 The non-linearity F satisfies Assumption IHJB1 with B be a singlton set 

Remark 4.5. Therefore, when the non-linearity F satisfies IHJB, we can drop the super script 
(5 and write F by 



where 



F(t,x,r,p,n/,i>) ■= inf x, r,p, 7) + l a (t,x,r,p,j,ip)} 



C a (t,x,r,p, 7 ) := ±Tr[(a a ) T ](t,x) 1 + b a (t,x)p + c a (t,x)r + k a (t,x), 



and 

l a (t,x,r,p,ip) := / (ip(x + T] a (t,x,z))-r-t { \ z \< 1} ri a (t,x,z)-p)v(dz). 
In this case, the non-linearity is a concave function of (r,p, 7, 

Assumption IHJB+ The non-linearity F satisfies IHJB2 and for any 5 > 0, there exists a 
finite set {ai} i= ^ 1 such that for any a E A: 



\k a - fc Q *|oc+yi(?f - Tf*)(;z)\l,di>(z) 



inf <^ \a a - a^loo + lfc" - 6 Ql |oo + |c a - c^U + l^ - £^|oo + /l(?? Q - rf*)(-, z)\i,dv(z) > <S. 

l<i<M s \ 



Remark 4.6. The Assumption IHJB+ is satisfied if A is a compact separable topological space 
and <7 Q (-), b a (-), and c a {-) are continuous maps from A to C fe 2 ' ([0, T] x ]R rf ); the space of bounded 
maps which are Lipschitz on x and ^-Holder on t and rj a (-) is continuous maps from A to |(/? : 

[0,T] x R d x R£-> r| j Ri Itp^z^uidz) < 00}. 

Theorem 4.2 (Rate of Convergence). Assume that the final condition g is bounded and Lipschitz- 
continuous. Then, there is a constant C > such that 
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" < C\hk + h6 2 K + he~ 3 + hU K + + h~i L zl< Jz\ 2 u{dz 



under Assumption IHJB, 

,K,h 

" I{\z\<k}\ 

• under Assumption IHJB+, 

-C (/l 1 / 10 + hT06 K + h^T K + fr& ! M < K} \z^v{dzj) < V - V K ' h 

In addition, if it is possible to find Kh such that 



lim k/j =0, limsupfo 4 ^^ < oo and limsup/i 2 / \z\ 2 v(dz) < oo, (4-14) 

h-+0 ft->0 J0<\z\<K h 



0<\z\<K h 

then, there is a constant C > such that 

• under Assumption IHJB, v — v Kh ' h < C/i 1 / 4 . 

• under Assumption IHJB+, — C/i 1 / 10 < v — v Kh ' h . 

Next example shows the case where the conditions of the above Theorem on the choice of k is 
satisfied. If it is not the case in some situations, it does mean that the rate of convergence is less 
than what is proposed by Theorem 14.21 

Example 4.1. For the Levy measue 

v{dz) = l K d\z\~ d ~ 1 dz, 

one can always find such that the condition of Theorem is satisfied. In the other words, it is 
always enough to choose Hh such that 

limsup h~2K h < oo. 

4.1. Convergence. We suppose the all the assumptions of Theorem 14 . 1 1 holds true throughout this 
subsection. 

We first manipulate the scheme to provide strict monotonicity by the similar idea as in Remark 
3.13 and Lemma 3.19 in |10j . Let u K,h be the solution of 

u K > h (T,-)=g and u K ' h (ti,x) = T Ki h[u K,h ](ti,x), (4.15) 

where 

T KA [^](t,x):=E^(t + h,X^ x ' K ^ + hF K>h {t,x,V h iP,iP(t + h,-)) (4.16) 



and 



F K>h (t, x, r, p, 7, VO =sup inf { -a ' 13 • 7 + b a >P ■ p + (c°^ + 6 K )r + e e ^ T ~^ k a ^(t, x) 

a P I 2 



+ 



{0fy l ^)-r-l M<l] n^{z)-p)v{dz) 
{M>«} 



Remark 4.7. Assumption Inf— Sup is also true if we replace 3% by 

JZ>P[ip](t, x) = l -a a ^ ■ D 2 cP(t, x) + b a >P ■ D(f>(t, x) + (c a 'P + 6 K )<t>{t, x) + e *»( T -*) k a >P (i, x) 

(Of ">\tj>(t,-))(x) -4>{t,x) - V a ^(t,x,z) ■ D<P(t,x))v{dz). 

{\z\>k} 

The proof is straight forward. 
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We have the following Lemma which shows that for proper choice of 6 K the scheme (|4,15p is 
strictly monotone. 



Lemma 4.1. Let 9 K be as in {4.&fy and ip and tf) : [0,T] 



be two bounded functions. 



Then: 



Proof. Let / := ip — ip > where ip and ip are as in the statement of the lemma. For simplicity, we 
drop the dependence on (£, x) when it is not necessary. By Assumption IHJB1 and Lemma (|3.1|) . 
we can write: 



- K.h I 



T #Cih fo>]=E[/(t + /i,X h )] 



+/i inf sup J?" [if>] (t + h,x)- inf sup J*' p [<p] (t + h,x)\, 
V a P a J 

where (f)(t,x) := K[<fi(t, X'£)] for cf> = ip or ip. Therefore, 

TK,hW\-T Kth [<p] > E[f(t + h,X h )] + h$%P$,ip](t + h,x,x), 
where is defined by 

^ [iP,<p](t,x):=M [</>](*,*)- sup J?? [ip] (t,x'). 

a p 

By Asumption Inf— Sup, there exists (a*,/3*) so that 

T^M-T^iv?] > E[f(t+h,x h ))+h(j*[^}{t+h,x)-j:mt+h,x)y 

Observe that by the linearity of Jk'^ ' , one can write: 

J^[4>](t + h,x) = E[j^[<p](t + h,X h ) . 
By the definition of Sk an d Lemma 2.1 in |lUj . 



T K M]-T K M > E 



1 „. 



f(X h ) (l + h(<g>f >m + B K + b?>fr.(v' 1 



h 



T,-lW h W^-hI d _! 



/?2 



-cr 



+w>r , '"* ,l (/) 



where 6^ = 6^ - / {1>w >^ a %)K^) and = - / {|z| > k} K^)- 
Therefore, by the same argument as in Lemma 3.12 in [10], one can write: 



T K ,ftW - T K ,^] > E 



/(*/*) 1 



i(6r^) T (a^^)-6r^))+W)f^' 1 



(/) 



where 



.4 



a*„S* 



(4.17) 
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Therefore, by positivity of / and Assumption IHJB1, one can deduce: 

T^l-T^M > HE f(X h )U + 6 K -j(bf>PY(a a *>P*)-bf>» 
By the choice of 6 K in (|4.8j) . we have 

T Kth [^)-T Kjh [cp] > 0. 
Then, sending e to zero provides the result. □ 

The following Corollary shows the monotonicity of scheme 14.31 
Corollary 4.1. Let tp,ifi : [0,T] x M. d — > R be two bounded functions. Then: 



In particular, if Kh satisfies then 

V < V> T Kh)h [<p] < T Khth [i/>] + ChEM -<p)(t + h, Xf > Kh )] 

for some constant C . 

Proof. Let 8 K be as in Lemma [4.1l and define <p K (t, x) := e 9K ^ T ~^ f(t, x) and ip K (t, x) := e dK ^ T ~ t ^(t, x) 
By Lemma |4.H 

T Kth [<p K ] < T KA [ip K }. 
By multiplying both sides by e~ dK ^ T ~ t \ we have 

( e -'»ft(l + 6 R h) - l)E[ V (t + + T^h^] 

< (e-t*h(i + 0^) _ i) E [i;(t + fc, If '")] + T Ki/l [# 

So, 

T« |h [^] < (e"^(l + - l) E[(V> - p)(t + h, Xj; x ' K )] + T^]. 
But, e"^ h (l + - 1 < -f^ e -»«\ So, 

T K ,,M < - e ^e-^ h E[^- V ){t + KX^ K )]+T K! M- 
which (I4.7P provides the result. □ 

In order to provide a uniform bound on v K ' h , we bound u K,h with respect to 6 K as in the following 
Lemma. 

Lemma 4.2. Let ip and ip : [0, T] x R rf — > R 6e too L°°— bounded functions. Then 

\T Kth [ip] - T K M}^ <\<p- V|oo(l + (C + 6 K )h) 
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where C = sup a g jc"'^^. In particular, if g is L°°— bounded, for a fixed k the family (u K,h (t, -))h 
defined in fi4.3\ ) is L°° — bounded, uniformly in h by 

Proof. Let / := <p — ip. Then, by Assumption Inf Sup and the same argument as in the proof of 
Lemma 14.11 



T Kth [tp] -T Kth [^]<E 



f(X h ) (l - a" 1 • a Q *^ + h(\A^* | 2 + + 9 K 

v{dz) - -(b a *^* - [ V a *' (3 \z)u{dz)) T (a a *'l 3 *y 

k\ 4\ Ju>\z\>k\ ' 



'{\z\>k} 

:[b a *'P 



V ^(z)u(dz))))+ht>f >\f) 
{1>\z\>k} J 



where A^ ^ is given by (|4.17p . On the other hand, 



< |/|oo / Hd 
J{\z\>k} 



Therefore , 

T K ,h[<P]-T*M < 1/looE 



I -or 1 - a a *^ + h(\Af>P*\ 2 + c a *>P* + 9 K 



Ub a *>F - [ v a *>P*(z)v(dz)) T (a a *^)-(b a *>P* - [ ri a *'P*(z)v{dz))) 

J{i>\z\>k} Jn>\z\>K\ 



{1>|2|>k} 



By Assumption IHJB1 and (fCgj) , 1 - a" 1 • a Q * ,r and 



rf* \z)v{dz)\ (a a *'P*)-(b a *^ 

{1>\z\>k} ' V 



rf*^{z)v[dz) 
{l>|z|>«} 



are positive. Therefore, one can write 



T K , h [ip] - T Kt h[ip] < |/U( 1 - a' 1 ■ a a *^ + h[E[\A% ' p ' \ z ] + c a ~^ + 9 K 



(4.18) 



1 

~ 4 

But, Notice that 



b a 



r, a *^{z)v{dz)) {a a *^)-(b a *t 

{1>|z|>k} ' ^ 



if' p \z)u{dz 
{i>l*l>«} 



E[\A?* ' r | 2 ] = h- 1 a~ 1 -a a *^ 



+ 



4V Jn>\z\>K\ ' v 



'{1>|z[>k} 



r] a ,p (z)u(dz 

{i>kl>*} 



By replacing E[\Af ,r | 2 ] into (^T5|l . one obtains 



T^M-T^M < 1/1^(1 + h(c a *'P* +0 K )) 

< \f\oo(i + (c + e K )h), 
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with C = sup a a |c Q '^|oo. By changing the role of ip and %p and implementing the same argument, 
one obtains 

\T K!h [<p]-%ML ^ \f\oo(i + (c+e K )h). 



To prove that the family (u K,h )h is bounded, we proceed by backward induction as in Lemma 
3.14 in [10] , By choosing in the first part of the proof tp = u K,h (ti+i, .) and ij) = 0, we see that 

\u K ' h (U,-)\oc < hCe "<- T -V + \u K > h (t i+1 ,-)\ 00 (l + (C + e K )h), 

where C := sup a ( g |fc Q '^|oo. It follows from the discrete Gronwall inequality that 

\u K > h (ti, .)|oo < (C(T - U) + \gUe {C+eK){T - h) . 

□ 

Define 

v K ' h := e- e ^ T -^u K ' h . (4.19) 
Next Corollary provides a bound for v K,h uniformly on k and h. 
Corollary 4.2. v K ' h is bounded uniformly on h and k, and 

\v K ' h - v K ' h \oo < K9^h for some constant K. 
If also, Kh satisfies fij. 7\ ), then 

lim \v Kh ' h -v Kh > h \oo = 0. 
Proof. By Lemma 14.21 for fixed k, we have: 

lu^Oloo < (C + lgUe^^. 

Therefore, 

\v"' h (t,.)\oo < (C+Hooje ^-*). 
For the next part, define u K,h (t,x) = e^ K ^ T ~^v K ' h (t, x). Direct calculations shows that 
-n,h = e e K h^ _ 0k q e r_«, fc ^ + ^ + h F K h ^ Xj v h u K ' h , u K ' h (t + h, ■) 

By an argument similar to Lemma 3.19 in |10| . we have 



|(u"' h -«"•*)(*, .)|oo < ^V'^ + Mloc (4.20) 

+ (1 + (C + 9 K )h)\(u K > h - u K > h )(t + h, OU, 
where C is as in Lemma 14.21 By repeating the proof of Lemma 14.21 for u K,h , one can conclude, 
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So, by multiplying 14.201 by e eK<yT ~ l \ we have 



2 ' 

+e- e * h (l + (C + K )h)\(v K > h - v K > h )(t + h, -)|oo, 



for some constant C. Because e ^(l + (C + 9 K )h) < e , one can deduce from discrete Gronwall 
inequality that 

\(v K ' h -v K ' h )(t,-)\oo < K9 2 K h, 
for some constant K independent of k which provides the second part of the theorem. □ 

We continue with the following consistency Lemma. 

Lemma 4.3. Let tp be a smooth function with the bounded derivatives. Then for all (t, x) E 
[0,T] x R d : 

hm - = -(£ <p + F(-,(p,D(p,D <p,<p(t, ■))) {t,x). 

(t' ,x')->-(t,x) 
(h,c)->(0,0) 
t' + h<T 

Proof. The proof is straightforward by Lebesgue dominated convergence Theorem. □ 



To complete the convergence argument, we need to proof the the approximate solution v Kfl,h 
converge to the final condition as 

Lemma 4.4. Let satisfy 7\ ), then v Kh)h is uniformly Lipschitz with respect to x. 

Proof. We report the following calculation in the one-dimensional case d = 1 in order to simplify 
the presentation. 

For fixed t G [0,T — h], we argue as in the proof of Lemma [4.21 to see that for x, x' £ M with x > x'\ 



u K ' h (t,x) -u K ' h (t,x') = E[(y, K >\t + h,X**)-u*\t + h,jt t ' x ') > )] 



+h inf sup J? ,P [u K ' h }(t + h,x)- inf sup J°^[u K < h ](t + h,x') 

\ a P a p 

u^ h {t + h, X^ x ) - u K ' h {t + h, X*'*'))} 



< E 

+h ( sup J?? [u^ h ] (t + h,x)- inf [u K M (t + h, x) ] . 

\p « ; 

Observe that by (14. 6j) . one can write 

u^fozj-u"^^') < E\(v«> h (t + h,jt t ' x ) -u K ' h {t + h,X^ x ')\ 

+. 



-h($ a >P{ U K> h ,u K > h }(t + h, x, x')) , 
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where <& is defined in the proof of Lemma 14.21 By Assumption Inf Sup, there exists (a*,/3*) such 
that 

^•^•[^^(t + ^x,^) = J?*'P*[^](t + h,x)-JZ*>FitM](t + h,a/). 
Therefore, 

u K ' h (t,x) -u K ' h (t,x') < ^u K ' h {t + h : X^ x )-u K ' h {t + h,X^ x ') 

+h (j* [W 1 ] (t + /i, x) - J K * [W 1 ] (i + h, x') ' 

For the other in equality we do the same except that when we 

u"' h (t,x)-u"' h (t,xf) < A + hB + hC, 

where 



A := E 



u K > h (t + h,X^ x ) - u K > n (t + h,X^ x '))] 



r t,X\ „.K,h 



+h ( jf'F [u^ h ](t + h,x)- Jf'f [u^ h ](t + h, x) 



with u K ' h (y) = u K ' h (y + x' - x), 



B := jf>F[u*fi](t + h,x)-jf>F[u*> h ](t + h,x?), 



and 



C := uf^'^u^t + h, -))(x) - uf^'^u^t + h, -))(x'). 
We continue the proof in the following steps. 



Step 1. 



C = h~ l E 



u R > h {t + h,X*> x )-u K > h {t + h,X 



where X*> x := x + Y%=i V°*'^ {x, Z t ) with Z iS are i.i.d. random variables distributed as — — - 



Step 2. By the definition of J K 



1 



B = -(a a *'P*(x)-a a *'P*(x'))V 2 h u K ' h (t + h,x') + (&f >?* (x) - V? {x'))V l h u K < h {t + h,x') 
+(c a *' /3 *(x) - c a *<P*(x'))Vlu K > h (t + h,x') + k a *'P* (x) - k a *>P\x'), 

where 6k ,/3 (x) := b a,/3 (x) — f^ 1> ^ z ^ >K y / q a ^(x,z)iy(dz). On the other hand, 



P£ = E 



Du^ h (t + h,X x h ')(^a-\x': 



k-i 



, for k = 1, 2. 
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So, 



B < E 



+(bf^(x)-bf^(x'))Du K ' h (t + h,X%') + (c a *>P*(x) - c a *^(x'))u K > h (t + h,X£ 
+f a *^\x) - f a *^(x'). 
Step 3. By the definition of , one can observe that 



jf'F[u*> h \{t + h,x)- Jf^[u^]{t + h,x) 



a*, /3* r~ K ,hl 



la<**>P*(x)SW+b* K (x)6(l)+c* K (x)5(0) 



where c* and 6* are defind in the proof of Lemma 14.11 and 



D k u K > h (t + h,X%)-D K u K ' h (t + h,X%) 



k„.K,h/ 



for k = 0,1,2. 



By Lemma 2.1 in |l(Jj . for k = 1 and 2 



E 



(>>^ + fr, X£) - u^\t + h, Xf)) H k h (t, x) 

a k (x) 



+u*> h (t + h,X%)H£(t,x)(l 



x' \ Tjk i 



a k (x') 



E 



u^{t + h, X x h ) - u K > n (t + h,X x h )) H k (t, x) 



+Du K ' h {t + h,X%)(^] a{x') ( a~ k (x) - a~ k (x r 



Therefore, one can write 



A < E 



u K > h {t + h,X%) -u K ' h {t + h,X£ 

x (l - a* + a*N 2 + frc* + b* K NVh"j (x) 
+hb* K {x')Du K ' h (t + /i, X%')a(x') {a-\x) - a^ix')) 
+a*{x')Du K ' h (t + h, xf)VhNa{x') {a~ 2 (x) - a~ 2 {x')) 



where a* := \a a *^\ a* := ^a' 1 ^*'^, c* := c a *^\ c* := c* + 6 K , and 6* := 6° 
Step 4. By dividing both sides by x — x' and taking the limit we have: 

Du K ' h {t,x) < E Du K ' h {t + h,Xl)(^l + hfi' K + y/ha'N + J K ^ 

x(l-a* + a*N 2 + hc* K + b* K NVh) 



«*,/?* 



+Du K > h (t + /i, X* h ' x ) (l+v*h+ f* h ) N%] + CeW-Qh, 
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where J Kjh := f {lzl>K} rj(z)J([0,h],dz), J* h :-- 
with intensity A K := Jn z \ >K \ v(dz). 
Let L t := \Du K ' h (t,-)\oo. Then 



Jn z t >K \ v'( z )J([Qi h], dz), and Nfi is a Poisson process 



E 



Du^ h {t + h,X; : x )(l + ^h + J* h )NZ < L t+h Ch(\ K + \':) 



where \'* := fn z \ >K \ rf*(z)v{dz). Let G := N + -^-\fh. By the change of measure 



:= exp 



4 2/ 



we have G ~ A/"(0, 1) under Q and one can write 



Du"> h (t,x) < ^^Du^\t + h,X^(l + h(Ji' K -^) + Vha , G + J^ 



„* to*) 2 



x 1 - a* + a*G^ + h(c* 



+h((b* K )' - b* K - -~-f) + {^^Ya- 1 - / 

+L t+h Ch (X K + X'*) + Ce 0K F-%, 



Step 4. Notice that 1 — a* + a*G 2 + /i(c* — ^p-) is positive and therefore, one can take m q, Z ] as 
a density for the new measure Q^. So, 



Du K ' h (t,x) < E l 



^ Du* fc (t + fc, X£) ( (l + h(p! K - ^f) + v^'G + J K)ft 



+L* +/l Cft (A* + A' K *) + Ce 9l < T -%. 



So, 



. a' bta 



+ Z-Hh (b* K y-b* K U -- U jf + i^)'*" 1 -^^ )VhG 



2 e qZ 



1 



1 + h(fi' K - ^) + v^a'G + J K , h 



a 2 / V2 
+Zt+ A C7i (A K + A',*) + Ce e « T -%. 



(7- 



Notice that 



E 



dP 



{Du^{t + h,Xt)f\ < L 2 t+he M\{b>?h). 
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On the other hand, 



1 + htf K - ^) + Vho'G + J K>h ) + Z- 1 (h((b* K )> ~bf-- h -^ 



+ 



cr 



E 



+ [^a a *^Ya- 1 -a a '^^)VhG] 

By calculation of the right hand side of the above equality, one can observe that all the terms of 
order \fh vanish and we have: 



A l + h{fi' R - ^) + Jho'G + J K>h ) + Z- l (h({bl)i ~bf-- h ^- 



(l(a^^)V- 1 



Kd' bio 



C 1 + h c* + 



* \2 



(K) 

4a* 



a z 



Therefore, by the choice of Kh, for h small enough we have 



L t < L t+h exp(^-h(C + 0K h - bK h *aa' + (b*K h 
< L t+h exp(h(C + 9n h ))+Ce e ^ T - t) h 



' - ^ - %^ + 2A Kfe + 2A'*^))+Ce^( T - i )/ l 



By discrete Gronwall inequality, 

Therefore by definition of v K ' h , we have 

\Dv^\ < eW-QQDgloo + CiT-t)). 



Lemma 4.5. Let Kh satisfies ( f^,7| j, th 



en 



limv K ' h (t,x) = g(x). 
Proof. We follow the same notations as in the proof of the previous Lemma and write 
u K ' h (t,x) = E\u«> h (t + h,jt t ' x )\ + hmfsupj^[^~ h ](t + h,x) 

L J a p 

< E\u K ' h (t + h,X^ x )\ + h sup J^[u^ h ](t + h,x). 



Observe that by (|4.6p . one can write 

u K ' h (t,x) < E\u"' h (t + h,X t > x )\ +h($ a >P[^~ h : ti\{t + h,x,x'j) + /tsup|/ 



□ 
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By Assumption Inf Sup, there exists (a*,/3*) so that 

u K ' h (t,x) < E\u K > h (t + h,X t > x )\ +hjf^[^ h ](t + h,x) + hC, 
where C := sup Q a |/ a '^|oo- Therefore, for any j = i,- ■ ■ , n — 1 one can write 



u* fc (t,,*2*) < u^{t 3+1 ,X%*){l-a* + a*G< + hC* 



ti .X 



+ hC. 



where a* := a*(tj, Xl*' x ), C* := (c* — )(tj,X^' x ) and G^s are independent standard Gauss- 
ian random variables under the new equivalent measure Q. By the consecutive use of the above 



inequality and the fact that 1 — a* + a^G 2 + hC* is positive, one can write 



n-l 



n-l 



u K > h (ti,x) < E Q g(X^ x ) 11 ( 1 - a* + a*G 2 + hC* ) + ChJ^ e 9 ' 



3 =i 



3=1 



Notice that in the above inequality we used the fact that 



1 - a* + a)G) + hC* 



i + hEY.[c*} <i + e K h. 



On the other hand, Z := U]=i (l - a j + a ) G ) + hC j) is positive there for fracZE®[Z] could be 
considered as a density of a new measure Q z with respect to P. Therefore, 

n-l 



u K ' h (ti,x) < E^[Z]E^ +ChJ3e M >. 



By the definition of one can write 



n-l 



+ e" 



e ^ . 



Therefore, 

^(ti,aO-<?(aO < e-^( T -^)E^[Z]EQ Z [| 5 (l^")- 5 (x)|] + C%(x)|(T - fc) + e^^C^T - U). 

Notice that g(x!f' x ) — g{x) converges to zero P-a.s. and therefore Q z a.s. as (U,h) — > (T, 0). So, 
by Lebesgue dominated convergence Theorem, 

limsup v K,h (ti,x) - g(x) < 0. 
(U,h)->(T,0) 

By the similar argument one can prove that: 



liminf v K,h (ti,x) - g(x) > 0, 
(U,h)^(T,0) 



which compelets the proof. 



□ 



Remark 4.8. By extending the above proof as in the Lemma 3.17 and Corollary 3.18 of [TO], one 
can proof that 

\v K > h (t,x) -g(x)\ < C(T-t)l 
Also, observe that by the similar argumnet as in |10| . v Kh,h is ^-Holder on t uniformly on h and x. 
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So, the approximate solution v Kh ' both satisfies the requirement of the convergence established 
in pQ and converges to a function v locally uniformly. Moreover, v is the unique viscosity solution 
of (|2.ip - (|2.2p . So, by Corollary 14.21 the same assertion is true for v K ' h . 

4.2. Rate of Convergence. For local PDEs, the rate of convergence of probabilistic numerical 
scheme relies on the approximation of the solution of PDE by regular sub and super-solutions, the 
consistency for the scheme and the comparison principle derived from strict monotonicity. One can 
approximate the solution of the local PDE by a regular sub-solution and a almost regular super- 
solution from up and down, respectively. These approximations are provided by a switching system 
and Krylov method of shaking coefficients. Next, we use the consistency Lemma 3.22 in [10] to 
produce inequalities for the regular approximations plugged into the scheme. Then by comparison 
principle; Proposition 3.20 in [10]; we obtain the bounds for the difference of approximate solution 
derived from scheme and regular approximate solution obtained from Krylov method and switching 
system. 

We continue this Subsection with establishing the same line of argument as in [10] for non-local 
case. The generalization of the method we used in [10] for the rate of convergence, is developed in 
[5] where the scheme needs to be consistent and satisfies comparison principle. Before, providing 
consistency and comparison principle result for the scheme (|4.15p . we show that truncation error 
could be handled by the Theorem of continuous dependence for (|2.1|) - (|2.2|) . More precisely, if v 
and v K are solutions of (|2.ip - (|2.2p and (|4.9p - (|4.10p . respectively; then by Theorem 5.1 in [6] 

\v-v K \oo < CI [ \z\ 2 u{dz) 

\J0<\z\<K 

Therefore, By choosing so that J 0< \ z \ <Kh \z\ 2 v(dz) < Ch? , one can just concentrate on the rate 
of convergence of v K,h to v K . 

We shift to v Kh,h which is is derived from the strictly monotone scheme (|4,15p and find the rate 
of convergence for v Kh,h . The following Corollary shows that this shift do not effect the rate of 
convergence. 

Corollary 4.3. For F which satisfies IHJB, F(t, x, 0, 0, 0, 0) = 0. Then, 

\v Kh ' h - v Kh ' h \ < ChQl . 

I I — Kfo 

In addition, if is such that 

- 2 

limsup/i^ < oo, 

then 

\v Kh ' h - v Kh > h \ < Chk 




Proof. The proof is straightforward by the proof of Lemma l4.2i 



□ 
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Form now on, we concentrate on the approximate solution v K,h which is obtained from strictly 
monotone scheme [415] through (|4.19j) . In order to provide the result, we need to use the consistency 
of the scheme for the regular approximate solutions. Then, the comparison principle for the scheme 
provides bounds over the difference between u K,h and regular approximate solutions. Let 

t? inn \ ^(t,x) - T Kih [iP](t,x) x 
TZ Kth [tjj\(t,x) := - h L ip{t,x) 

+F K (;^,D^,D 2 4>,iP(t,-))(t,x). 
Lemma 4.6. For a family {^ £ }o<e<i of smooth functions satisfying 

df°L>V < Ce 1 - 2 ^-^ 1 f or any ^ G N x N d \ {0}, 

where := Yli=i Pi> an< ^ C > is some constant, we have: 



C I he' 3 + h6 K E~ l + h^ K + e~ l I \z\ 2 u{dz)) , 



(4.21) 



|^,ftbdloo ^ R ( h > £ ) 
for some constant C > independent of k. If in addition 

limsup/i^ < oo and limsup\//i / \z\ 2 v{dz) < oo, 

/i->0 /i->0 J{\ Z \< K \ 

we have: 

\KK h ,h[<Pe]L ^ R ^ £ ) : = C (he^ + Vhe- 1 ). 
Proof. 1Z K ^\(p £ ] is bounded by 

sup^E^^t + h,Xl x ' K ) - <p £ (t,x)) + ^Tr [a a (D 2 Ve (t + h,Xj; x ' K ) - D 2 <p E (t,x)) 
+b a {Dip e (t + h, Xj; x ' K ) - D<p e (t, x)) + (e K + c a ){y £ {t + h, Xj; x ' K ) - <p e (t, x)) 

+I a [<p e ](t,x)-l% h [<p e }(t + h,x) 

For the Levy integral term by Lemma 13.31 we have: 

\l a [v £ ](t,x)-l« h [<p £ ](t + h,x)\ < CHDip^Vh + h f \z\v(dz)) 

V J{\z\>k} 

+h\d t D 2 ip £ \ 00 + \D\ e \ 

J{\z\<k} J 

< C (he- 3 + /i v / ^ + £~ 1 / \z\ 2 v{dz)\. 

V h\A<K} J 

By the same argument as Lemma 3.22 in [TU] all the other terms are bounded by he~ 3 except 
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which is bounded by O^he . The second assertion of the Lemma is straightforward. □ 

Next we need to have maximum principle for scheme 14.151 Note that Lemma 3.21 in pJJj holds 
true for scheme [4TT5] with A = 8 K + C and (3 > 9 K + C where C = sup Q \ c a \. Therefore, Proposition 
3.20 in [10] holds true for non-local case. More precisely, we have the following Proposition. 

Proposition 4.2. Let Assumption F holds true, and consider two arbitrary bounded functions if 
and ij) satisfying: 

h~ l (ip - ThM) < 91 and h~ l (if, - T h ty>]) > g 2 
for some bounded functions g\ and g 2 - Then, for every i = 0, ■ ■ ■ , n: 

(<P - Mfax) < e^ +c ^\(<p - ^)+(T, Oloo + (T- h)e^ +c ^ T -^\{ gi - 52 )+| 00 
where C\ > sup a \ c a \. 

The approximation of the solution of non-local PDE by the Krylov method of shaking coefficients 
and switching system is developed in [5]. [5] provides the result of rate of convergence of general 
monotone schemes for the non-local PDEs satisfying Assumption IHJB. However, they referred 
the regularity of the approximate solutions to the result of [6] where the approximate solution 
obtained from switching system could only be locally ^-Holder continuous on t. But, in the case 
of scheme (|4.3p . we need the solution of (|2.ip - (|2.2p be uniformly ^-Holder continuous on t. It is 
because we need the regular approximate solutions obtained from Krylov method and switching 
solution to satisfy (|4.2ip . Therefore, in the present work we need to rebuild Lemma 5.3 in [6] under 
the Assumption IHJB to obtain global ^-Holder continuous on t for the solution of the switching 
system. 

Therefore, we continue this subsection by introducing the switching system of non-local PDEs 
with the regularity result needed for the solution of this system. 
Let k be a non-negative constant. Suppose the following system of PDEs: 

max {— C x vi(t, x) — Fi (t, x, Vi(t, x), Dvi(t, x), D 2 Vi(t, x),Vi(t, •)) , Vi — M. l v} = 

«iCT, ■)=#(•), (4-22) 

where i = 1, • • • ,M and 

Fi(t,x,r,p,j,^) := inf {£ a (t,x,r,p, 7,7) +l a (t,x,r, p, 7,^)} 

aeAi 

£ a (t,x,r,p,7, 7 ) := ^Tr [a a (t, x) 7 ] + b a (t, x) ■ p + c a (t, x)r + k a (t, x) 
l a {t,x,r,p,j,^) := / (4>(t,x + ri a (t,x,z)) - r - t { \ z \< l} r] a (t,x, z) ■ p) dv{z) 
A4 l r := mmrj + k. 

We would like to emphasize that gis need to satisfy gi — A4 l g < where g = (51, • • • , <?m)- If each 
g-i = g then we obviously have gi — Ai l g < 0. 



30 



A PROBABILISTIC SCHEME FOR FULLY NON-LINEAR NON-LOCAL PARABOLIC PDES 



Existence and comparison principle result for the above switching system is provided in Propo- 
sition 6.1 [5]. Also, it is known from Theorem 6.3 in [5], that if (v 1 , ■ ■ ■ ,v M ) and v be respectively 
the solutions of (j4T22j) and (j2J]) - (fZ2|) with A = Uf^A and As are disjoint sets, then 



< v % - v < Cks for i = 1, ■ ■ ■ , M. (4.23) 

The regularity result for (I4.22P is provided in [6]. There, it is proved that (v l ) is Lipschitz with 
respect to x and locally l/2-H61der continuous with respect to t. For the proof of Theorem 14.21 
(V) should be uniformly l/2-H61der continuous with respect to t. The following Lemma provide 
the uniform l/2-H61der continuity for (v l ). 

Lemma 4.7. Assume HJB holds for each i and let (V) be the viscosity solution of \4--22 ). Then 
there exist a constant C such that for any i = 1, • • • ,M: 

Mi ^ c - 

Proof. Lipschitz continuity with respect to x is done in Lemma 5.2 in [6]. To obtain uniform 
1/2— Holder continuity with respect to t, the proof of Lemma 5.3 in [6] should be modified by using 
assumption HJB. 

Fix y 6 M. d , t and t' where t < tl '. For each i = 1, • • • , M, define: 



L 



i/>i(t,x) := A— e A(t ~ l) \x - y\ 2 + B(t'-t) + K(t' — i) + A — + v l {t' : y) 



Where L = and A, a and 7 will be defined later. Then: 

d t ^(t,x) = -\^(Ae A V-V\x-y\ 2 + B)-K 



So, 



2 

Dij>i(t,x) = \Le A{ t'-V{x-y) 
D 2 ^(t,x) = ALe^'-^irfxrf. 



-d t ^i + C a {t, x, ipijDipijD 2 ^) + X a {t, x, Vi, D^i) 
= -\L(Ae A V-V\x-y\ 2 + B)-K 



+\Le A ^'-^Ti [a a (t, x)] + 2\Le A{ ±' ~% a (t, x) ■ (x - y) + c a (t, x)if>i + k a (t, x) 
+ \Le A ^^ [ (\x + if(t, x, z) - y\ 2 -\x- y\ 2 - 21 {W < 1} r ? a (i, x, z) ■ (x - y)) dv(z). 



2 

By HJB, we choose K and A so that, 

KU < K, |6 a U < K, I^U < K,\k a \oo < K, K^ 1 < X < K 

\v\oo < K,\r] a (t,x,z)\ <K{lA\z\). 

Without loss of generality and with the similar argument as in Remark 3.7, we can suppose that 
for any a, c a < 0. So by choosing positive large D, there exists positive constants C\, C2 and C3 
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such that: 

-dtifn + C a {t, x, fa Difji, D 2 ^i) + l a (t, x, faDfa) 
< -\LeW-*) K ^-^Ix-yf-^x-yl+^B- C7 2 ) + C 3 . 

Therefore, choice of large B and D makes the right hand side negative. 

-d t tp i + C a (t,x,i; l ,Di; l ,D 2 i; i )+l a (t,x,i; i ,D^ i ) < 0. 

On the other hand, 

i>{t', X ) = ^(\\ X -y\2 + \-±) +V i( t ',y). 

Minimizing with respect to A, 

i/j(t', x) > L\x -y\+ v^t', y) > v\t', x). 
We can conclude that ipi is a super solution of (|4.22|) . So, by comparison Theorem in [5], 

Mt,v)>v%v)- 

So, 

L(XB{t' - t ) + X- 1 )+v i (t',y)>v i (t,y). 

Therefore, 

v i (t,y)-v\t',y)<CVt F ^t. 
The other inequality can be done similarly by choosing: 
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4>i(t,x) := -A 



L 



A(t'-t)\ x - y |2 _ Bit 1 - t)] - K{t' -t)-\- 1 - + v l (t', y). 



□ 



Remark 4.9. Notice that all the result of switching system is correct for (|2.ip - (|2.2p satisfying F 
by simply setting M = 1 and k = 0. 

Therefore, by [5] there are regular functions w£ and w 1 ^ which are respectively the regular sub- 
and super-solution of 

-C x u K (t, x) - F K (t, x, u K (t, x),Du K (t, x), D 2 u K (t, x),u K (t, •)) = 0, on [0, T) x R d , 
u K (T, •) = g, on G R d . 

where 

F K (t,x,r,p,~f,7p) := mf\C a ^(t,x,r,p,-f)+l^(t,x,r,p,'y,ijj)\ 
(one can replace sup inf by inf sup) where 
£ Q ^(t,x,r,p, 7 ) := ^Tr 



a a >Pa a >P T (t, s) 7 + b a '^(t, x)p + (c a '^(t, x) 



+ 
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and 



l^(t, x, r,p, 7, ip) := / 0(z + ^(t, x, z)) - r - l { \ z \< 1} v a ' P (t, x, z) ■ p)v{dz). 

J{\z\>k\ 



Therefore, by Proposition 6.2 and Theorem 6.3 of [5], Lemma 14.61 and Proposition 14.2 

(u K -u K ' h )(t,x) < (u K -y£ + u£ -u K ' h )(t,x) 



< Ce {e M) {T-t) ( + h£ -s + Mk£ -i + h ^/r K + £ -i [ \ z \2 

\ ^{|z|<k} 



v(dz) 



and 

(u K > h -u K )(t,x) < (u K ' h l -w K £ +w R £ -u K )(t,x) 



< cjoMw-t) L\ + h£ s + h0K£ -i + h ^- K +£ -i r | 2 



z\ v{dz) 
} 



Notice that v K (t,x) = e~ ^ T ~^u K {t,x). So, 

V n_-K,h < C[e + hE-* + M K e- s - +hy/0 K + e- L I \z\ A v(dz) 



and 



< C le + he- 3 + h9 K e- 1 + h y ^ + e~ 1 [ \z\ 2 
\ J{\z\<k\ 

v K > h -v K < C I +he~ 3 + h9 K e~ 1 + h v ^ + e^ 1 [ \z\ 2 u{dz)\ 
\ Al*l<«} / 



z\ 2 v{dz) . 

K] J 

On the other hand, because of (|4.14p and by Lemma (|4.2p . the second part of Theorem 14.21 is 
provided after choice of optimal e. 

5. Conclusion 

The algorithm is the first probabilistic numerical method for fully non-linear nonlocal parabolic 
problems. As in local case (|10j). it converges to the viscosity solution of the problem. A rate of 
convergence is known for the convex (concave) non-linearities. Also with the same argument as 
in Section 4 in [10] . Monte Carlo approximations of expectations inside the scheme do not affect 
the asymptotic results if enough number of samples would be used. The error analysis for MCQ 
shows that the appropriate approximation of jump-diffusion process with compound Poisson process 
could be applied in discretization procedure. The theoretical result is followed by some numerical 
examples which confirms the convergence of the scheme. 

On the other hand there are some features where the scheme is not implementable in non-local 
case, e.g. when the non-linearity is of HJB type. This could be the challenge of future works. The 
other open issue is to relax some assumptions. For example, relaxing the assumption of uniform 
ellipticity may be a future work. 
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